Potts model with q = 3 and 4 states on directed Small- World network 



(N 

o 

(N 
oo 



P. R. 0. da Silva, 1 F.W.S. Lima/'E and R. N. Costa Filho 2 'Q 

1 Dietrich Stauffer Computational Physics Lab, Departamento de Fisica, 
Universidade Federal do Piaui, 64049-550, Teresina - PI, Brazil 
2 Departamento de Fisica, Universidade Federal do Ceard, 
Caixa Postal 6030, Campus do Pici, 60455-160 Fortaleza, Ceard, Brazil 

Monte Carlo simulations are performed to study the two-dimensional Potts models with q = 3 and 
4 states on directed Small- World network. The disordered system is simulated applying the Heat 
bath Monte Carlo update algorithm. A first-order and second-order phase transition is found for 
q — 3 depending on the rewiring probability p, but for q = 4 the system presents only a first-order 
phase transition for any value p . This critical behavior is different from the Potts model on a 
square lattice, where the second-order phase transition is present for q < 4 and a first-order phase 
transition is present for q > 4. 
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It was conjectured by Harris [l[ that the sign of 
the critical exponent of the specific heat a determines 
whether spin systems are affected or not by randomness. 
For positive values of a the system with randomness 
or impurities has a critical behavior different from the 
pure system case. For negative values of a, on the other 
hand, the critical behavior of the system should be the 
same for both pure and impure cases. In particular, for 
two-dimensional regular lattices, the ferromagnetic Potts 
model with ostates displays first order phase transitions 
for q > 4 while the pure ferromagnetic three-state 

Potts model has a — 1/3, hence, according to the above- 
mentioned criterion we expect to find a different behavior 
for a random interaction system. However, Picco [4J] and 
Lima et al. studied this model with different type 

of disorder and did not find any relevant difference from 
the pure case. 

The q-state Potts model has been studied in scale-free 
networks by Igloi and Turban Q and depending on the 
value of q and of the degree-exponent 7 first- and second- 
order phase transitions were found. This model was also 
studied by Lima [l(| on directed Barabasi- Albert (B A) 
networks, where only first-order phase transition has 
been obtained for any (/-values with connectivity z = 2 
and z = 7 of the directed BA network. Here, we stud- 
ied the Potts model with q = 3 and 4 states. We also 
calculate the critical exponents ratio j3/v and 7/1/ for 
second-order phase transitions that appears due to the 
SW disorder. 

We consider the ferromagnetic Potts model with q = 3 
and 9 = 4, on directed small- world networks where every 
site of a directed small-world network of size N — L x L 
have spin variables a taking values 1,2,3 and 1,2,3,4 
for q = 3 and 4, respectively. With L being the side of 
a square lattice. In this network, created by Sanchez ct 
al. [ll[ (see Fig. 1), we start from a two-dimensional 
square lattice consisting of sites linked to their four near- 
est neighbors by both outgoing and incomplete links. 
Then, with probability p, we reconnect nearest-neighbors 




FIG. 1. Sketch of a directed small- world networks constructed 
from a square regular lattice in d — 2. Figure gently yielded 
by Juan M. Lopez from Sanchez et al. [111 ]. 



outgoing links to a different site chosen randomly. After 
repeating this process for every link, we are left with a 
network with a density p of SW directed links. Therefore, 
with this procedure every site will have exactly four out- 
going links and different (random) number of incoming 
links. The time evolution of this system is given by a 
single spin-flip like dynamics with a probability pi : 



1 



[l+exp(2£ 2 /fc s T)]- 



(1) 



The Hamiltonian of a ^-states ferromagnetic Potts 
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model can be written as 



H=-J J2 S*< 



(2) 



<hj> 



where S is the Kronecker delta function, and the sum 
runs over all neighbors of <7j. 

The simulations have been performed applying the 
HeatBath update algorithm on different lattice sizes: 
N = 64, 256, 1024, 4096, and 16384. For each system 
size quenched averages over the connectivity disorder are 
approximated by averaging over R = 40 independent re- 
alizations. For each simulation we have started with a 
uniform configuration of spins (the results are indepen- 
dent of the initial configuration). We ran 4 x 10 4 Monte 
Carlo steps (MCS) per spin with 2 x 10 4 configurations 
discarded for thermalization. 

In studying the critical behavior of the model using 
the HeatBath algorithm we define the variable e = E/N, 
where E is the energy of system, and the magnetisation 
of system M — (q. maxfni] — N)/(q — 1) , where rii < 
N denote the number of spins with 'orientation" i = 
1, q. From the fluctuations of e measurements we can 
compute: the average of e, the specific heat C and the 
fourth-order cumulant of e, 



u(T) = [< E >} av /N, 



C(T) = K 2 N[< e 2 > - < e > 2 } av , 



B(T) 



< e 1 > 
3 < e 2 > 5 



(3) 



(4) 



(5) 



the temperature can be defined as T = J/ksK, where ks 
is the Boltzmann constant. Similarly, we can derive from 
the magnetization measurements the average magnetiza- 
tion (m = M/N), the susceptibility, and the magnetic 
cumulants, 



m(T) = [< |m| >]„ 



X (T) = KN[< m 2 > - < \m\ >% 



Ui(T) = 



1 - 



< rrr > 



3 < m > 



(6) 



(7) 



(8) 



where in all the above equations < ... > stands for a 
thermodynamic average and [...] QU for an average over 
the 40 realizations. 

To verify the transition order for this model, we ap- 
ply finite-size scaling (FSS) (l2{ . Initially we search for 
the minima of the fourth-order parameter of Eq. (5). 
This quantity gives a qualitative as well as a quantita- 
tive description of the order of the transition 13j. It 
is known [141 ] that this parameter takes a minima value 



-Bmin at effective transition temperature T C (N). One 
can show [lj| that for a second-order transition limAr^^ 
(2/3 — .B m in) = 0, even at T c , while at a first-order transi- 
tion the same limit measuring the same quantity is small 
and (2/3 - B min ) ? 0. 

A more quantitative analysis can be carried out 
through the FSS of the C fluctuation C max , the suscep- 
tibility maxima Xmax and the minima of the Binder pa- 
rameter -Bmin- 

If the hypothesis of a first-order phase transition is 
correct, we should then expect, for large systems sizes, 
an asymptotic FSS behavior of the form 16|- 18| , 



a c + b c N + 



Xmax = a x + b x N + 



a B + b B /N - 



(9) 



(10) 
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if the hypothesis of a second-order phase transition is 
correct, we should then expect, for large systems sizes, 
an asymptotic FSS behavior of the form 



C = C reg + L a /»f c (x)[l + ...], 
m = L-e/»f m (x)[l + ...}, 
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L 1/V fu(x)[l + ...], 



(12) 



(13) 



(14) 



(15) 



where C reg is a regular background term, v, a, j3, and 
7 are the usual critical exponents, and fi(x) are FSS 
functions with 



(16) 



being the scaling variable, and the brackets [1 + ...] indi- 
cate corretions-to-scaling terms. Therefore, from the size 
dependence of M and \ we obtain the exponents (3/v 
and "f/v, respectively. The maxima value of susceptibil- 
ity also scales as L 1 ^ . 

For each value of q, we apply the finite size scaling tech- 
nique [12j . and the same procedure is done for systems 
with different number of sites N = 64, 256, 1024, 4096, 
and 16384. The critical temperature for infinite size sys- 
tem is estimated by using the fourth-order magnetization 
(Binder) cumulant. 

In Fig. 2, we show the dependence of the energy u and 
magnetization m on the temperature T, obtained from 
simulations on directed SW network with lattice size L = 
8, 16, 32, 64, and 128 and the rewiring probability p = 0.0, 
0.1, and 0.9. The shape of m(T) and energy u curve, for 
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FIG. 2. Display of the energy (top panel) and magnetisation 
(bottom panel) against temperature T for p — 0.0 (circle), 
p = 0.1 (square), and p — 0.9 (diamond). Here L=128 and 
we are considering the case when q — 3. 
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FIG. 3. Plot B e (T) versus T for: a) p = 0.1 and b) p = 0.9 
for different size lattices L — 8 solid line, L = 16 dotted line, 
L — 32 dashed line, L = 64 long dashed line, and L = 128 
dotted-dashed line. In all cases we used q = 3. 



the particular parameters used (N = 16384 and q = 3), 
suggests the existence of a second-order phase transition 
in the system for p = 0.0 and p = 0.1, and a first-order 
phase transition in the system for p = 0.9. The phase 
transition occurs at the value of the critical parameter 
T c . The energetic Binder cumulant as a function of the 
reduced temperature T is shown in Fig. 3 for p — 0.1 and 
0.9 and different lattice sizes (L = 8 to 128). From the 
figure one can see a typical second-order phase transition 
(for a large system B e (T) — > 2/3 ) and a first-order phase 
transition is observed for p = 0.1 and 0.9, respectively. 

In Fig. 4, the difference 2/3 — B m i n is shown as a 
function of the parameter 1/y/N for p — 0.1 and p = 0.9. 
For p = 0.1, a second-order transition takes place since 
the limAr^oo (2/3 — Bi !m i n ) — 0, even at T c . However, 
for p = 0.9 a first-order transition is observed, because 
one has (2/3 — B i]min ) ^ 0. 

We display the scalings for natural logarithm for the 
dependence of the magnetization m on inflection point 
at K = T C (L) and p = 0.1 for q = 3 in the Figure 5. 
The slopes of curves correspond to the exponent ratio 
P/u according to Eq. 13. The obtained exponents are 
fi/v = 0.24(5). The exponents ratio 7/1/ are obtained 
from the slopes of the straight lines with 7/1/ = 1.5(1) 
for SW, as presented in Fig. 6 and obtained from Eq. 
14. The results present a reliable indication in favor of 
the Harris criterium, error bars are only statistical, and 
much larger systems might give different exponents, also, 
the exponents ratio /3/v and 7/1/ obey the hyper-scaling 
law j/i/ + 2P/i/ =d. 
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FIG. 4. Plot of 2/3 - B min at T c as a function of 1/VN for 
several values of the system size N = 64 to 16, 384 sites for 
p = 0.1 (circles) and p — 0.9(squares). 



Next, we study the case where q — 4. In Fig. 7, as in 
the Fig. 2, we show the dependence of the magnetization 
to and energy u on the temperature T, obtained from 
simulations on directed with lattice size L = 8, 16, 32, 64, 
and 128 with (LxL = N) sites and the rewiring probabil- 
ity p = 0.0, p = 0.1, andp = 0.9. The shape of m(T) and 
energy u curve, for a given value of N = 16384 sites and 
q = 4, suggests the presents of the second-order phase 
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FIG. 5. Display of the magnetisation at the inflection point 
versus the size system L p = 0.1, q = 3. 
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FIG. 6. Logarithmic plots of the suscepbility at T c versus the 
size system L for p — 0.1 , q — 3. 
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FIG. 7. The same plot of Fig. 2, but now for q = 4. 
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FIG. 8. The same plot of Fig. 5, but now for q — 4. 



transition in the system for p = 0.0, but also suggests 
the presents of the first-order phase transition in the sys- 
tem for p = 0.1 and 0.9. 

In Fig. 8, as in the Fig. 5, we plot the difference 
2/3 — B m i n as a function of the parameter 1/yN for 
different probabilities p — 0.1 and p — 0.9 . Unlike the 
q = 3 case, for both values p = 0.1 and 0.9 a first-order 
transition is observed, because (2/3 — Bj, m i n ) 7^ 0. 

In conclusion, we have presented simulations for Potts 
model with q — 3, and 4 states on directed SW network. 
The disordered system is simulated applying the Heat- 
Bath Monte Carlo update algorithm. The Potts model 
with q — 3 does display a second-order phase for rewiring 
probability p — 0.1, with exponent ratio j3/v — 0.24(5) 



and 7/1/ = 1.5(1) that are different of the Potts model 
on a regular lattice, where, the specific-heat exponent 
a = 2/3 is a good candidate for a change of the critical 
exponents, that agree with the Harris criterium [l[ and 
obey the hyper-scaling law 7/1/ + 2[3Jv = d and for case 
of p = 0.9 we have a first-order phase transition. In the 
case g = 4 both values here studied rewiring probability 
p = 0.1 and 0.9 present a first-order phase transition as 
showed in the Fig. 7 and 8, that again agree with Harris 
criteria. In summary, the behavior of Potts model for 
q = 3 and 4, here studied, is due to the directed links of 
the SW networks, where can have short and long range 
interaction. 
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